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ABSTRACT 

Using small-angle X-ray scattering, we determined 
the three-dimensional packing architecture of the 
minichromosome confined within the SV40 virus. In 
solution, the minichromosome, composed of closed 
circular dsDNA complexed in nucleosomes, was 
shown to be structurally similar to cellular chromatin. 
In contrast, we find a unique organization of the 
nanometrically encapsidated chromatin, whereby 
minichromosomal density is somewhat higher at 
the center of the capsid and decreases towards the 
walls. This organization is in excellent agreement 
with a coarse-grained computer model, accounting 
for tethered nucleosomal interactions under viral 
capsid confinement. With analogy to confined liquid 
crystals, but contrary to the solenoid structure of 
cellular chromatin, our simulations indicate that the 
nucleosomes within the capsid lack orientational 
order. Nucleosomes in the layer adjacent to the 
capsid wall, however, align with the boundary, 
thereby inducing a 'molten droplet' state of the chro- 
matin. These findings indicate that nucleosomal 
interactions suffice to predict the genome organiza- 
tion in polyomavirus capsids and underscore the 
adaptable nature of the eukaryotic chromatin archi- 
tecture to nanoscale confinement. 

INTRODUCTION 

The formation of viruses from proteins and nucleic acids is 
a striking example of complex self assembly. Many viruses 
form icosahedral, almost spherical capsids, inside which 
nucleic acid is packaged. Bacteriophages, for example, 
package linear double-stranded DNA (dsDNA) using 



motor proteins, thereby achieving dense, highly ordered 
and pressurized DNA organization (1-7). Other viruses 
package linear single-stranded RNA (ssDNA) through 
direct electrostatic interactions between RNA and capsid 
protein that nucleate capsid assembly (8-10). These 
viruses typically show high density of RNA lining the 
interior of the viral capsid (11-14). 

The genome of simian virus 40 (SV40), a non-enveloped 
icosahedral virus and member of the polyomaviridae (15), 
poses additional constraints to packaging. The viral 
genome is a 5.2-kb circular dsDNA, wrapped around 
histone octamers, forming a chromatin-like structure 
with ~20 nucleosomes, termed the 'minichromosome' 
(16-18). Unlike cellular chromatin or minichromosomes 
isolated during the infection cycle (18), minichromosomes 
packaged within the mature virion lack histone HI 
proteins, most likely facilitating condensation. The viral 
capsid surrounding the viral chromatin comprises three 
viral-encoded proteins VP1, VP2 and VP3. In total, 
72 VP1 pentamers arrange in a T = 1 icosahedral lattice 
with ~50-nm external diameter and ~36-nm internal 
diameter (19,20). Long, flexible carboxy-terminal arms 
tie the pentamers together in three different types of inter- 
actions, allowing gaps between the subunits. 

During assembly, capsomers are added around the 
minichromosome (21) leading to its compaction. In this 
process, the minichromosome must overcome the frustration 
arising from internucleosomal (predominantly steric) inter- 
actions, as well as confinement by the capsid walls. The 
three-dimentional, higher-order structure of the closed 
circular minichromosme within the capsid of polyomaviruses 
(as well as papillomaviruses) confined inside this nanometric- 
sized cavity remains unknown. In contrast, cellular chroma- 
tin is known to be present as ~10-nm fiber in 'beads on 
string' or 'open', transcriptionally active structure, or alter- 
natively as ~30nm 'closed', often inactive fiber configur- 
ation. Although the in vivo prevalence and relevance of the 
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30-nm chromatin fiber is still under debate (22,23), a number 
of different higher-order structures for this configuration 
have been suggested by several authors, including solenoid, 
zigzag, superbead, and other architectures (24-30). The un- 
certainty in structure appears to be related to the particular 
nature of the nucleosome arrays studied, and to the sensitiv- 
ity of the 30-nm fiber to changes in external conditions, wit- 
nessed both experimentally and by high-resolution and 
coarse-grained simulations [see ref (30) for a recent extensive 
review]. Interestingly, once released from the virion, the 
SV40 minichromosome assumes chromatin structures rem- 
iniscent of cellular chromatin, which depends on pH (18) 
and salt concentration (17). 

An early theoretical model of the packaged 
minichromosome based on volume restrictions (31) 
concluded that the chromatin may fit into the capsid 
with little configurational freedom. While at that time 
the precise shape and size of the nucleosome and capsid 
volume were ambiguous, the model highlighted the strong 
nucleosomal confinement within the SV40 capsid. 

Dubochet and co-workers attempted to capture the 
packaged minichromosome structure by rapidly vitrifying 
samples following lysis of virion particles (32). However, 
once the minichromosome was released from the confines 
of the virions, the minichromosome quickly assumed a 
different structure. Freshly prepared minichromosomes 
(containing 20-25 nucleosomes) at close to physiological 
ionic strength (130mM NaCl) were found condensed into 
globules, 30 nm in size or more. At low ionic strength 
(<10mM NaCl), the minichromosomes opened, first 
into lOnm filaments and then into nucleosome-strings. 
The globule-like structure was explained by a liquid 
drop model, whereby nucleosomes assemble through 
stacking interaction to form distorted (solenoid, 30 nm) 
or fragmented 10-nm filaments. These structures cannot 
represent the tight sphere-like confinement of the 
packaged minichromosome, yet are consistent with the 
interactions in solution of nucleosomal core particles 
(NCPs), as demonstrated by detailed experiments (33,34). 

NCPs result from digestion of inter-nucleosomal linker 
DNA and can be approximated as discs of 5nm height 
and 10 nm diameter (35). Depending on their concentration 
and on solution parameters, such as pH and osmotic 
pressure (exerted by PEG polymer) (33,34), NCPs show 
complex phase behavior in solution, and can assemble 
into well-defined macroscopic liquid crystalline order (36). 
For example, similar to phases found for discotic or oblate 
ellipsoids that are favored by steric packing considerations 
(37), NCPs at densities found inside the virus, and at 
physiological conditions, form columnar or even higher- 
order phases (36,38). Under physiological salt concentration 
(~150mM), NCPs show net attraction, as determined by 
osmotic pressure measurements that yield the so-called se- 
cond virial coefficient (33). This attraction may further 
facilitate the formation of condensed columnar phases (36). 

To gain insight into the interplay between minichro- 
mosome confinement and molecular interactions, we used 
a combination of solution small-angle X-ray scattering 
(SAXS) experiments and coarse-grained computational 
modeling. Our results suggest that, when averaged over a 
large ensemble of viruses, the SV40 minichromosome within 



the virion adopts configurations consisting of somewhat 
higher chromatin density in the center and lower mean 
densities for larger radii. Good agreement between simula- 
tions and experiments allows us to rationalize the forces and 
interactions leading to the observed organization, and to 
predict structural features that have not yet been resolved 
by experiments. Importantly, we find that weak, non- 
specific and steric interactions are sufficient to account for 
the minichromosomal structure within the SV40 capsid. 

MATERIALS AND METHODS 

Preparation of wt SV40, empty capsid and DNA-VLPs 

SV40 stocks were prepared as previously described (39) 
and separated into full (DNA-containing) and empty par- 
ticles as detailed in SI sec. la. Virus-like particles (VLPs), 
containing dsDNA (DNA-VLPs), were prepared as previ- 
ously described (40,41) and detailed in SI sec. lb. 

SAXS 

The methodology followed for measuring the diffraction 
of viral particles has been previously described (42,43) and 
is detailed in SI sec.lc. 

To model the SV40 samples, a form factor of multiple 
concentric spherical shells with smoothly varying electron 
density (ED) profiles, represented by hyperbolic tangent 
functions: 

Ap(r) = 0.5 J Api + [( A /°'+i " A <°<-)] tanh[j ( -(r - R,)] J 

(1) 

was fit to the data, as described (42,43). Ap(r) indicates the 
radial ED contrast with respect to the solvent (buffer). 
The index i represents the rth layer in the sphere, with 
an outer radius, R h an ED contrast, Ap„ and linked to 
the subsequent layer (i+ 1) by a slope s t . Ap N = 0 and 
corresponds to the solvent ED contrast. To numerically 
solve arbitrary ED profiles, the form factor of the smooth 
ED profiles, Ap(r) were calculated by discretizing the 
profiles and transforming them to a sum of Heaviside 
step functions that can be analytically solved (42,43). 

An improved fit (T = 1 model) to the scattering data 
from wtSV40 was obtained by replacing the protein 
capsid layer (the third layer in Table 1) with spheres of 
smoothly varying ED contrast (Equation 1 with N = 2), 
located at the center of each pentamer, represented by: 
n f N-l 

/=0 I i=l 

x tanh[si(|r| - i?,)] J+0.5{ Apcapsid - Ap Capsi d (2) 

72 

x tanh^copridd?! - «)]} ® ^8(7-/ x R„) 

n=\ 

The first term in the expression represents the ED profile 
inside the capsid as a sum of multiple concentric spherical 
shells as described in Equation 1. The second term 



Nucleic Acids Research, 2013, Vol. 41, No. 3 1571 



Table 1. Parameters for the best fit models of the SAXS data (presented in Figure la-c) for wt SV40, empty capsids and VLP (DNA-VLPs) 



Model parameter 


wt SV40 


Empty capsid 


DNA-VLPs 


b Ri [nm] 


9.2 


16.0 


16.7 


u Ap, (pi) (e/nm 3 ) 


40 (373) 


13 (346) 


2 (335) 


sj (nm ') 


0.2 


0.3 


0.4 


h R 2 (nm) 


18.1 


24.4 


24.9 


a Ap 2 (p 2 ) (e/nm 3 ) 


37 (370) 


43 (376) 


55 (388) 


s 2 (nm ) 


0.7 


7.2 


1.6 


b i? 3 (nm) 


24.4 






a Ap 3 (p 3 ) (e/nm 3 ) 


49 (382) 






sj (nm ) 


0.5 






Background 


0.0739 ? ~ 2 - 43 + 0.0343 


0.076 ? - 2 - 84 + 0.566<T 0 - 36 


0.417?~ 196 + 0.031 



"'The error in the ED is ±2. 
b The error in the radii is ±0.1 nm. 



represents the arrangement of the 72 VP1 pentamers in the 
capsid shell. This term results from a convolution between 
a sphere, with an outer radius a and ED Ap Capsid , and^ a 
sum of delta functions shifted by a set of vectors fxR„. 
The R„ vectors represent the locations of the center of 
mass of the 72 VP1 pentamers derived from the crystal 
structure of wtSV40 (20). The fitting parameter / enables 
us to change the radius of the particle. 



Monte carlo simulations 

Using the Metropolis algorithm (44), we performed 
Monte Carlo simulations in the canonical (NVT) 
ensemble, whereby nucleosomes were modeled as oblate 
ellipsoids confined within a spherical shell, as 
parameterized by experiments and further detailed in SI 
sec.ld-g, Supplementary Figures S1-S4. 

The attractive potential between the nucleosomes is 
known to be determined by the salt concentration in 
solution. In our simulations, we assume that the 1:1 salt 
concentration in solution is set by the parameter X that, 
more generally, represents the Lagrange multiplier 
associated with constraints set by solution conditions. 
This parameter is used in the standard metropolis algo- 
rithm to determine transition probabilities as 



P(a -+ b) 



1, E b < E a 

e -{E b -Ea)/Xk B T £ > £ 



(3) 



Where P(a — > b) is the transition probability at any 
step from a state with energy E a to a state with 
energy E b . 

Variation in X in the simulations leads to a modulation 
of effective forces experienced between interacting par- 
ticles, through its impact on the second virial coefficient, 
B 2 , as determined in the simulations (Supplementary 
Figure S5, sec. 2a). Variations in experimentally 
determined B 2 with salt concentrations are available for 
NCPs in solution under various concentrations of salt 
(33). By mapping experimental B 2 with salt concentration 
and simulated B 2 with X, allows us to derive X versus salt 
concentration. 



RESULTS AND DISCUSSION 

SAXS measurements show that SV40 has higher 
chromatin density in the capsid center 

Because SV40 is nearly spherical, as seen by the crystal 
structure (20), we fit the radially integrated X-ray scatter- 
ing intensities as a function of the magnitude of the 
momentum transfer vector, q, (Figure 1) with a model of 
concentric spherical shells with a smoothly varying radial 
ED profile (Materials and Methods Equation 1). This ap- 
proximation should yield a good fit at lower q values (cor- 
responding to the larger elements of the icosahedral 
capsid), whereas at higher q values (~0.6nm _1 , corres- 
ponding to ~10.5nm in real space), where molecular 
details become important, the model is not expected to 
fit as well. Note, that SAXS measures the averaged struc- 
ture of many virions, rather than the structure of a single 
specific virion. 

As in earlier scattering studies (45^17), in addition to 
subtraction of the buffer intensity, we subtracted, from the 
buffer-subtracted data, a series of empirical power laws 
that pass through the minima of the scattering oscillating 
intensities and coincide with the highest available </-data 
(q > 3.5 nm -1 ), as previously described (42). This add- 
itional background originates in concentration fluctu- 
ations of our multi-component mixture (full particles, 
partially assembled particles, capsid proteins and nucleic 
acids) that lead to a (/-dependent SAXS. Concentration 
fluctuations are Lorentzian-like (peaking close to q = 0), 
hence contribute to power-law (q~ 2 ) background. In 
addition, reactants (DNA and/or capsid proteins) and 
partially assembled structures may contribute to the back- 
ground. The underlying assumption of our background 
subtraction is that most of the oscillating scattering 
signal comes from fully formed structures with uniform 
dimensions (i.e. the distribution of dimensions, e.g. 
particle radius, throughout the sample is very narrow). 
Our previous studies showed that both wt SV40 and 
DNA-VLPs are monodispersed [see transmission 
electron microscopy micrographs shown in Figure 6B 
and C in (40)]. We have therefore based the modeling on 
that premise. Moreover, had the sample been 
polydispersed, for example, by less than a nanometer in 
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wt SV40 7=7 model 




0.2 0.4 0.6 0.8 1 0 

q [nm" 1 ] 
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DNA-VLPs 



0.2 0.4 0.6 0.8 1.0 

q [nm' 1 ] 




0.2 0.4 0.6 0.8 1.0 

q [nm 1 ] 




Figure 1. Radially integrated solution SAXS intensities (symbols) versus the magnitude of the momentum transfer vector q of (a) wt SV40, (b) empty 
capsid (c) VLPs containing dsDNA (DNA-VLPs) (d) wt SV40 T = 7 model. For (a-c), the solid red curves are the best fitted form-factor models of 
multiple spherical shells with a smoothly varying radial ED profiles, represented by hyperbolic tangent functions (Equation 1) and shown as insets to 
each plot. For (d), the solid red curve is the best fitted J ' = 7 form-factor model (Equation 2). The average radial ED is shown as insets. The dashed 
blue curves show the assumed power-law background functions (see text and Tables 1 and 2). The best fitted parameters of the three models (a-c) are 
presented in Table 1 and for T = 7 form-factor model (d) in Table 2. The R 2 values for the models of wt SV40, empty capsid, DNA-VLPs and wt 
SV40 with the T = 7 model are 0.997, 0.979, 0.999, 0.996, respectively, (e) Representation of the fitted geometrical T = 7 model for wtSV40 
(Equation 2). The red core represents the ED of the confined minichromosome, and the blue spheres represent the 72 VP1 pentamers building 
the virus capsid. The transparent spheres enable to view the capsid interior. Doted spheres mark the pentamer vertices of the icosahedron. The blue, 
yellow, red and green bars correspond to inter-pentamer distances of 9.0, 10.3, 9.7 and 8.5 nm, respectively, as the best fitted parameters establish. 
The SAXS data were measured at the ID02 beamline at the ESRF synchrotron (Grenoble). The experimental q range is q — 0.07556-4.423 14 nm~'. 



the radius of the particle, the deep minima seen in Figure 1 
would have been filled by a significant amount, as 
demonstrated in (43). 

The best model fit for the wt SV40 virus (Equation 1) 
indicates high ED at the external spherical shell and in the 
center (Figure la). This finding suggests a high ED 
associated with the proteinaceous outer capsid layer 
(with typical ED values of 382 ± 10e/nm 3 , Table 1) and 
a uniform radial distribution of the nucleosomes within 
the capsid. In the model, the SV40 capsid's outer radius 
is 24.4 ± 0.1 nm, and the inner radius is 18.1 ± 0.1 nm 
compared with 24.7 nm and 17.9 nm determined for the 
crystal structure, respectively (20). These data allow a 
precise (sub nm) size estimate of the SV40, under physio- 
logical solution conditions. Furthermore, our results agree 
well with the density profiles determined in cryo-EM 
measurements of wt SV40 (48) (Figure 2a). Unlike the 
present study, the earlier crystal structure did not 
provide any information about the internal organization 
of the minichromosomes. 



As expected, the spherical shell model fit does not 
account for features at higher q values, reflecting the 
icosahedral capsid deviations from a sphere. The dis- 
crepancy is seen beginning at g>0.65nm _1 (Figure 
la). The deviation of the double maximum (at 
(7«0.7nm -1 ) arise from the arrangement of the cap- 
sid proteins, which, based on our more detailed T = 1 
model (Equation 2 and Figure Id), leads to the double 
hump observed in the SAXS data, corresponding to a 
mean distance between spheres of 9.5 nm (Figure le). 
This value is the typical distance between pentamer 
centers. In SV40, the pentamer cores, 6-7 nm in size, 
are separated by gaps due to their unique binding via 
long carboxy-terminal arms (20). The capsid's outer and 
inner radius, as measured by the detailed T = 7 model 
(Table 2), is similar to the radius calculated with the 
spherical shell model (Equation 1 Table 2 and Figure 
la). Moreover, the detailed T= 1 model and the spher- 
ical shell model both indicate a uniform distribution of 
the nucleosomes within the capsid with a possible 



Nucleic Acids Research, 2013, Vol. 41, No. 3 1573 
Table 2. The best fitted values for the scattering form factor of wtSV40, based on the T = 7 geometrical model (Equation 2 and Figure Id) 



Capsid parameters Value Inner core parameters: Value 

Parameters for the capsid spheres Parameters of the inner spherical shells 



a ED (e/nm 3 ) 391 b Layer 1- outer radius (nm) 8.3 

b Radius (nm) 4.7 b Layer 2- outer radius (nm) 17.1 

Slope (nnT 1 ) 3.6 Slope layerl (mrT 1 ) 0.5 

Icosahedron grid radius (nm) c 21.2 Slope layer2 (nm -1 ) 1.2 

Inner radius of capsid (nm) d 16.5 a E.D. layer 1 (e/nm 3 ) 369 

" l E.D. layer 2 (e/nm 3 ) 368 

Total radius of inner core (nm) 17.1 



Background function 0.0739 x 9" 243 +0.0 3 43 



The capsid parameters describe the size, position and ED of the spheres building the T = 7 icosahedral capsid (second term in Equation 2). The inner 
core parameters describe the averaged ED profile of the minichromosomal structure (the capsid interior given by the first term in Equation 2). 
An empirical power law background function was added to the form factor model (dashed blue carve in Figure Id), as discussed in the text. 
"The error in the ED is ±2(e/nm 3 ). 
b The error in the radii is ±0.1 nm. 

The icosahedron radius is the distance between the center of the virus and the center of each sphere that is placed at each of the 12 vertices of the 
T — 7 structure. 

d The inner radius is equal to the difference between the icosahedron grid radius and the radius of the spheres representing the VP1 pentamers. 



slightly higher (statistical presence) nucleosome density 
in the center of the capsid. 

Empty capsids, separated from wt SV40 by equilibrium 
sedimentation in a CsCl gradient (SI sec Id), generate 
scattering data that is best fit by low ED within the 
capsid (Figure lb and Supplementary Figure S9a). The 
difference between the ED of the outer shell and the 
center is significantly larger for the empty capsids than 
for wt SV40. This result substantiates that the viral 
minichromosome accounts for the high ED at the center 
of the wt SV40 particles. Yet, the ED of the inner volume 
of the empty capsid is higher than the solvent ED, possibly 
because of unidentified molecules of the host cell trapped 
within the empty capsids. These can also account for the 
thicker wall of the empty capsid (8.4 nm) compared with 
the wall thickness of the wt SV40 (6.3-6.6 nm, for the 
T = 7 model the wall width is calculated at half 
maximum). 

In the absence of nucleosomes, dsDNA packs close to 
the capsid shell 

Synthetic VLPs containing dsDNA (DNA-VLPs) are 
assembled in vitro from recombinant SV40 capsid VP1 
proteins (produced in insect cells) and closed-circular 
(histone-less) dsDNA (40) (SI sec. lb). These DNA-VLPs 
do not contain VP2 and VP3, which are present in both wt 
SV40 and empty particles purified from CV1 primate cells. 
The scattering curve of these DNA-VLPs shows that their 
wall ED is higher than that of wt SV40, whereas its thick- 
ness resembles that of the empty capsid. Their ED at the 
center is lower than that of both empty capsids and wt 
SV40 (Figure lc and Supplementary Figure S9b). We in- 
terpret this finding as follows: owing to large persistence 
length of dsDNA (~50nm), nucleosome-free DNA within 
the DNA-VLPs is forced against the outer wall, thus 
adding to its thickness (compared with the wt virus) and 
leaving the center clear. This explanation is supported by 
our earlier findings indicating that during assembly in vitro 
of DNA-VLPs, the plasmid DNA must undergo compac- 
tion (40). This DNA organization is reminiscent of that 



found in dsDNA bacteriophages, where the DNA concen- 
trates close to the viral capsid owing to its large persist- 
ence length (1,49). In contrast, in wt SV40, the 
nucleosomes condense the DNA, relieving the force that 
drives DNA to the capsid wall and allowing it to spread 
across the interior of the virion. These results indicate that 
the higher-order nucleoprotein structure of the minichro- 
mosome dictates nucleosome distribution in the capsid. 

Coarse-grained simulations describe viral chromatin 
packaging 

To investigate the possible structure of the 
minichromosome within the SV40 capsid, we modeled 
the confined chromatin arrangement using a coarse 
grained model, and validated the model by comparing 
with our SAXS experiments. Different levels of coarse 
graining have previously been used to model various 
aspects of oligonucleosomal arrangement (24,29,50-59). 
These simulations were used to investigate the role of 
histone tails, salt concentration, nucleosome repeat 
length, entry-exit angle of DNA at the histone surface 
and more. Specifically, simulations indicate that the 
longer the DNA linker between nucleosomes, the larger 
the diameter of the resulting chromosome fiber (56,57). 
Moreover, experiments (27) have shown that for the 
calculated average linker length of the SV40 minichro- 
mosome, ~90bp (see SI sec.le), the expected fiber 
diameter is ~45nm, raising the question of how the 
minichromosome can fit into a ~36nm diameter capsid. 
Additional simulations indicate that at higher salt concen- 
tration (200 mM), the attractive electrostatic interactions 
between the nucleosomes are more dominant, thus leading 
to more condensed structures (58). None of these models, 
however, investigated the chromatin arrangement within 
confined cavities, similar to the conditions existing within 
the SV40 virion. 

Here, we use a simple mesoscopic model, wherein viral 
chromatin is composed of oblate ellipsoids, each repre- 
senting a nucleosome. Based on experimental data (35), 
the ellipsoids have a ratio of short to long radius of 
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Of/a e — 0.5, and interact through Gay-Berne-type 
potentials (60,61). This orientationally dependent inter- 
action potential accounts for excluded volume repul- 
sions and longer ranged attractive forces between 
nucleosomes; these forces were parameterized based on 
experimentally determined interactions between NCPs in 
solution (33,34) (SI sec. Id). In our model, ellipsoids 
(histone octamers) are tethered via harmonic potentials, 
representing linker DNA stretches between nucleosomes. 
In the model, we define the average distance between the 
nucleosomes, r 0 = 30 nm (or 3ct with <7 = er e =10nm) 
where 30 nm corresponds to ~90bp (SI sec.le), and 
their average sliding distance along the DNA as ± 1 5 nm 
(SI sec.le). 

The SV40 icosahedral capsid is modeled as a spherical- 
bounding shell interacting with nucleosomes via short 
range (steric) repulsive forces and long range, mostly elec- 
trostatic, attractive forces, acting normal to the shell 
surface (SI sec. If). Guided by parameters derived from 
our SAXS experiments (Table 1), we set the effective 
internal radius of the capsid sphere (the maximum 
allowed radius for the minichromosome) at 20 nm. 
Based on previous reports, we assumed N = 20 ellipsoids 
in a capsid (16-18). 

Molecular configurations are generated using the 
standard Metropolis algorithm (44). Thus, our simulations 
predict the average nucleosomal distribution of an 
ensemble of many virions at their free energy minima, but 
cannot determine the chromatin order of a specific virion or 
the kinetics of nucleosome movement within an individual 
capsid. The effect of salt concentration in the bulk is rep- 
resented through the parameter X ('Materials and Methods' 
section). Values of X = 1.5-1.75 correspond to 100-150- 
mM bulk salt concentration (the physiological and our ex- 
perimental ionic strength), as calibrated from NCP solution 
data, and considering the reported error in those measure- 
ments (33) (see also SI sec 2d). Higher or lower X values 
correspond to lower and higher salt concentrations, respect- 
ively. Even though we associate X with salt concentration, it 
can similarly be correlated through similar mappings with 
other thermodynamic variables such as temperature, pH or 
osmotic stress. We assume room temperature in all simula- 
tions, as in our experiments. 

To gain insight into the forces and constraints that de- 
termine the nucleosomal order observed in simulations, 
our strategy was to consider several reduced models. In 
each of these models, a subset of interactions was either 
modified or completely turned off, allowing us to evaluate 
its impact on nucleosomal arrangement by comparison 
with the 'full model' (FM) described in the earlier 3 
paragraphs, which considered all interactions. Thus, to 
follow the impact of nucleosome connectivity, we used a 
reduced model that modifies the DNA-tethering forces 
between nucleosomes (tethering model, TM), and a 
model that assumes that the nucleosomes are free to 
slide along the DNA (no-tethering, NT). Similarly, to 
estimate the interactions of nucleosomes with the 
boundary, we studied a model that, in addition to free 
sliding, included a purely repulsive force with the capsid 
walls (NT and no-attractive-boundary, NT-NAB). 



Experiments and simulations show similar chromatin 
densities in the virion 

In excellent agreement with the SAXS experiments, the 
FM indicates high average nucleosome densities in the 
center of the capsid that decreases towards the periphery 
(compare FM and SAXS, Figure 2a). In comparing 
simulations with experimental ED profiles, we present 
the simulated average local ED distribution along the 
radius of the capsid d(r), normalized such that 
J cap d(r)Aitr 2 dr = 1. The remarkable agreement between 
experiments and simulations suggests that the packaging 
constraints of the nucleosomes, together with simple inter- 
actions considered in the simulations, are the main driving 
forces leading to the structure and order of chromatin in 
the SV40 virus. At salt concentrations of 50-300 mM, cor- 
responding to A > 1.5, the change in radial density is 
minimal (Supplementary Figure S8, X = 1.5, 1.75 and 3). 
In contrast, at lower X (corresponding to salt concentra- 
tion >>300mM), the nucleosomes become more ordered 
and start assuming a solid-like organization 
(Supplementary Figure S8, X = 0.7 and 0.3). We 
conclude that modest deviations from physiological con- 
ditions are predicted to be insufficient to change the or- 
ganization of the packaged minichromosome. 

Rigid nucleosome tethering would lead to strong 
nucleosome stacking at the capsid wall and less 
accumulation in the capsid center 

Modifying tethering interactions let us estimate the extent 
of packing frustration that nucleosomal connectivity 
ensues. We compared three variations of tethering inter- 
actions for nucleosomal connectivity (Figure 2b). In the 
first tethering model (TM-1 nm), nucleosomes are tethered 
by rigid harmonic linkers with an equilibrium distance of 
30 nm and a potential well-width corresponding to mean 
fluctuations (sliding distance) of ± 1 nm at thermal energy, 
k B T. This model corresponds to rigid DNA linkers with 
the known DNA persistence length (~50nm at physio- 
logical conditions) and with no nucleosome sliding 
allowed. In the second model, the nucleosomes are 
tethered strongly, with mean nucleosome-nucleosome dis- 
tances set at 30 nm and mean fluctuations of ±10nm 
(TM-lOnm). This corresponds to the minimum sliding 
indicated experimentally for the few SV40 nucleosomes 
experimentally found to be strongly localized (62), (see 
SI sec.le for further details). In the third model, nucleo- 
somes are freer to move (slide) with respect to one 
another, with the average distances (30 nm) but fluctu- 
ations of ±35nm (TM-35nm). For comparison, we also 
present results for untethered nucleosomes (NT). 

In contrast to the experiments (Figure 2a, SAXS), simu- 
lations of the strongly tethered nucleosomes (Figure 2b, 
TM-1 nm and TM-lOnm) show particles tightly packed 
against the capsid internal wall in a dense (hypertrophic) 
layer. This finding suggests that this model incorrectly 
captures nucleosomal mobility, in agreement with 
previous experimental determination of nucleosome 
phasing in the virus (62). Those experiments showed that 
strongly localized nucleosomes are widely separated from 
one another, with long stretches of loosely localized ones. 
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Figure 2. Chromatin local density d(r) in simulations as would be 
observed along the radius of the capsid. (a) For three models at 
X = 1.75: Blue curve, FM; Green, in the absence of tethering (NT); 
Red in the absence of tethering and lacking attractive forces to the 
boundary (NT-NAB). Also shown are the best model fit for the 
SAXS data for SV40 (taken from the inset to Figure Id), normalized 
to the FM density at low r (orange dashed line); the spherically 
averaged density profile of SV40 measured by Cryo-EM adapted 
from Shen et al. (48) and normalized to FM density at low r 
(magenta dashed line), (b) d(r) at X = 1.75 for green, NT model; 
Black, strong tethering well (minimal sliding, TM-1 nm); Red, weaker 
tethering (less sliding, TM-lOnm); Navy, weak tethering (maximal 
sliding, TM-35nm); Blue, the FM. See text for additional details. 

Indeed, SAXS experiments fit much better with simula- 
tions of more loosely tethered nucleosomes that assume 
nucleosome sliding of ± 15nm (compare FM and TM-10 
Figure 2b with SAXS, Figure 2a). In contrast, when 
tethering is weak or absent (Figure 2b, TM-35 nm or 
NT, respectively) we find a higher nucleosome density in 
the center of the capsid and a reduced density near the 
capsid walls with respect to our experimental data. Even 
though the structure of the minichromosome (closed 
circle) is probably more constrained than linear chroma- 
tin, simulation results for circular and linear chromatin, as 
expected, are identical (within statistical error). To fully 
appreciate the possible importance of the topology of 
circular chromatin, additional forces absent from our 
model, such as supercoiling and the resulting torsion 
(63), should be considered. However, our current stage 
of understanding of chromatin structure does not allow 
to reliably incorporate the topology. Nevertheless, the FM 
well reproduces the available experimental results, from 
both SAXS and cryo-EM, supporting its use as a 
minimal model for the SV40 minichromosome. 
Moreover, previous modeling studies suggest that the 
overall higher-order structure of chromatin is dominated 
by inter-nucleosome interactions, and is rather insensitive 
to details of the linker free energy (64). 

Attractive capsid walls increase nucleosome density at the 
hypertrophic layer, while maintaining the higher density in 
the capsid center 

The capsid protein VP1 possesses a DNA-binding domain 
(65). Moreover, the calculated electrostatic potential on the 
capsid interior shows 'salt-and-pepper' speckles of posi- 
tively and negatively charged residues (Supplementary 
Figure S3), which could lead to attractive forces with the 
minichromosome. Supplementary Figure S3 shows an 



internal view of the capsid crystallographic structure (left) 
and the corresponding electrostatic potential as calculated 
using Poisson-Boltzmann theory (right) (66). To evaluate 
the role of these attractive forces between the capsid and 
the minichromosome, we simulated untethered nucleo- 
somes that lack an attractive force with the boundary 
(NT-NAB), as described in SI sec. If. These conditions 
result in high nucleosomal density in the capsid center, 
even higher than those found when an attractive 
boundary is used, (compare NT-NAB with NT Figure 
2a). The lack of attractive forces with the boundary also 
results in lower local density in the hypertrophic layer than 
when attractions are included (Figure 2a). Clearly, even in 
the presence of moderate attractive forces with the wall, the 
density maximum at the center of the capsid persists, re- 
flecting dominance of the inter-nucleosomal attraction and 
spatial confinement (Figure 2a, NT). 

An additional reduction of the FM substitutes Gay- 
Berne potentials (ellipsoids) with nucleosomes modeled as 
spheres, interacting through Lennard-Jones potentials, and 
trapped within a spherical repelling boundary (SI sec. 2b). 
As in the FM, this reduced model indicates strong particle 
density in the center of the cavity (Lennard-Jones Model- 
LJM in Supplementary Figure S6). The fact that both el- 
lipsoids and spheres lead to elevated local density in the 
cavity center supports our conclusion that this particle ar- 
rangement, found in SV40, reflects dominance of the 
confining shell wall. Notably, this confinement effect is 
abolished when the system is confined in only one dimen- 
sion, such as a simple liquid between two walls, where 
density in the center is never the global maximum. 

Contrary to nucleosomes in bulk solution, nucleosomes 
within SV40 are disordered 

The good correspondence of simulations with the experi- 
ments (Figure 2a, FM and SAXS) provides validation to 
the simulations, indicating that the major forces respon- 
sible for particle ordering in the virion are appropriately 
accounted for. With this, we turn to explore additional 
information that is hardly accessible by available experi- 
mental techniques. Specifically, we use the model to 
predict possible nucleosomal orientations within the 
capsid. The reference point for our analysis is the well- 
studied state of nucleosomes in the bulk. 

To make links to bulk experimental measurements of 
NCPs in solution (38), we simulate untethered ellipsoids 
under periodic boundary conditions (PBC) representing 
the bulk. To follow the possible ordering of particles, we 
use the orientational order parameter, 



where 8 is the angle between the ellipsoid director and the 
director defining the average orientation of all particles, z 
(as detailed in SI sec. 2c). Although the simulated systems 
considered here are finite (TV = 100), we borrow termin- 
ology from the descriptions of orientational phases and 
their transitions in infinite systems, as follows. At high a 
(low salt concentration), where interactions between par- 
ticles is weak, the ellipsoids are in an isotropic-like regime 
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and do not show preferred alignment (low S). Yet, as X 
decreases and attractive interparticle interactions become 
dominant, we find an orientational transition, and the el- 
lipsoids preferentially align along a common axis, corres- 
ponding to a bulk nematic phase (nematic-like regime). 
Concomitantly, 5 increases and particles enjoy stronger 
favorable contacts (asterisks mark transitions points in 
Figure 3a, PBC). Similarly, as X is further lowered (cor- 
responding to very high salt concentration), the system 
becomes increasingly more ordered and undergoes an add- 
itional transition into a columnar-like regime and finally 
to a solid-like regime (compare snapshots in Figure 3b). 
These results are consistent with simulation results of 
oblate ellipsoids in the bulk (67). 

Reported studies of unconfined NCPs in aqueous solu- 
tions (38) indicate that nucleosomes order in a columnar 
phase at densities and ionic strengths similar to those of 
the native virus (~150mM salt and ~200mg/ml NCPs 
corresponding to X = 1.5-1.75 and nucleosome reduced 
density p Nuc m 0.6, respectively; for definition of pn uc , 
see SI sec. 2a). As expected, S values indicate that for 
our bulk simulations at X = 1.5-1.75 untethered ellipsoids 
are ordered in the nematic-like regime, or situated close to 
the isotropic-nematic transition (Figure 3a, PBC N = 100, 
shaded area). In contrast, for N = 20 nucleosomes within 
a spherical shell (representing SV40) at the same A's, the 
nucleosomes assume a 'molten droplet' state, whereby 
particles are on average orientationally disordered (FM, 
NT-NAB, NT Figure 3a shaded area). This effect is found 
regardless of DNA tethering or boundary attraction, 
indicating that confinement is sufficient for reducing the 
order of chromatin. 

The small number of nucleosomes is sufficient to change 
the salt concentration threshold for the phase transition 

We find that the order-disorder transition (isotropic-like to 
nematic-like) occurs at a specific 'critical' value X c . 
Following several numbers of particles at similar densities 
indicates that, the lower the number of particles, the lower 
the determined X c (Figure 3a and Supplementary Figure 
S7). This is analogous to the known finite-size effect, 
which manifests as a shift in the critical phase transition 
temperature of other known order-disorder transitions 
(68-70). Here, these shifts are the result of the limited 
system size that restricts the correlation length for orienta- 
tional fluctuations, which must necessarily diverge to allow 
a second-order transition (71). We conclude that, in accord- 
ance with this finite-size effect, the smaller number of nu- 
cleosomes is sufficient to decrease (compared with bulk 
solutions) the value of X c required for the order-disorder 
transition. This finding suggests that at physiological ionic 
strength, the confined nucleosomes are in a disordered 
state, whereas NCPs in the bulk (N >> 100) are ordered; 
for example, for N = 300 we find X c x 1.9. 

Capsid confinement can further decrease nucleosomal 
ordering 

Although for the bulk, there are two clear orientational 
transitions (market by asterisk, Figure 3a, PBC), for the 
confined, untethered ellipsoids, we find only one clear 



transition from isotropic-like to nematic-like organization. 
As X is further lowered (corresponding to higher salt), the 
particles under confinement become increasingly more 
ordered, moving to a columnar-like arrangement regime 
(Figure 3a, snapshots), but not to the solid-like regime 
(Figure 3b, left). Moreover, while for the ellipsoids 
simulated in the bulk at low X, we find that S approaches 
1 (Figure 3a, PBC) for the ellipsoids confined within a 
spherical shell, the highest order parameter reached, re- 
gardless of DNA tethering or boundary attraction, is 
5 = 0.8 (Figure 3a FM, NT and NT-NAB). This effect 
is also found for all simulations of ellipsoids confined 
within a spherical boundary, at various numbers of par- 
ticles, densities and cavity diameter (Supplementary 
Figure S7). Thus, we conclude that the confining spherical 
boundary is sufficient to induce disorder in the system and 
hinder the solid-like regime. 

Nucleosomes in the hypertrophic layer align with the 
capsid walls 

To resolve the boundary effects on nucleosomal order, we 
define the radial order parameter, E, in terms of the ellips- 
oid alignment with the spherical shell's radius vector (SI 
sec. 2d). From the variation of E with distance from the 
boundary, we find that, for both low and high X, the el- 
lipsoids within the 'wetting layer' near the boundary show 
an average radial order parameter close to E = 1 
(Figure 4, red and black curves). This finding suggests 
that particles mostly lie flat against the boundary, thus 
reducing their excluded volume interactions with the 
bounding wall. However, only at low X, the order param- 
eter 5 for ellipsoids outside the wetting layer (defined as 
r< 1.5a) indicates that particles preferentially align with 
respect to one another significantly (compare red and 
black horizontal line, Figure 4). Thus, particle ordering 
is affected by the balance between inter-particle inter- 
actions favoring high S and wall-particle interactions 
favoring high E, as demonstrated in snapshot of 
Supplementary Figure S7b. These competing forces 
become more prominent at low X, where both 5 and E 
are higher. As a result, the system can exist in the 
ordered columnar-like regime (S = 0.8) arrangements, 
but owing to the requirement to also order with the 
boundary, cannot reach the high degree of order seen in 
the solid-like regime arrangements (S = 1). 

The balance of wall and bulk interactions determine 
nucleosomal order and the critical point 

Simulations show that X c decreases when interactions with 
the boundary become more dominant, specifically, by 
increasing the nucleosome local density in the hyper- 
trophic layer as a results of an attractive potential at the 
boundary, or of DNA tethering, (compare FM, NT with 
NT-NAB, Figure 3a). This decrease in X c indicates that 
orientational frustration with the wall becomes stronger 
in comparison with inter-particle attractions, even at 
lower X. Therefore, the model predicts that a higher 
ionic strength (corresponding to screening of repulsive 
interactions) is required before the confined ellipsoids 
can overcome the aligning tendencies with the boundary, 
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Figure 3. Nucleosomal order in simulations, (a) Order parameter S 
versus X (associated with the concentration of salt) at a constant 
density of p Nuc « 0.6. Navy is for PBC N = 100; Black, PBC N = 20; 
Red, NT-NAB model; Green, NT model; Blue, FM. Top right insets 
show snapshots of the system from two different perspectives for the 
FM at X - 0.25, taken from simulations conditions marked in plots by 
blue square. The highlighted area marks the range 1.5 < X < 1.75, corres- 
ponding to the range of physiological salt concentrations. Dashed lines 
mark the average S for 100 or 20 randomly oriented ellipsoids. Black and 
orange asterisks mark the critical point for two consecutive transitions 
(isotropic-like to nematic-like and nematic-like to columnar-like, respect- 
ively), calculated via the derivative of S. (b) Snapshots of 100 ellipsoids 
under bulk-like conditions (i.e. under PBC) at constant density, 
Pnuc ~ 0.6. Left, snapshot of the system for X = 0.1, taken from simula- 
tions conditions marked in plots by navy square. Right, snapshot of the 
system for X = 1 with additional snapshot of a slice through the center in 
the at plane. All snapshots are oriented so that the eigenvector of the 
highest order parameter is defined as f . 




Figure 4. Nucleosome orientational order along the capsid radius from 
simulations. Full dotted lines follow the radial order parameter E(r) 
versus distance from the boundary, for the FM. Short straight lines 
indicate the average order parameter (S) for all the ellipsoids that are 
located away from the boundary (/• < 1.5a-), shown for X = 1.75 (black) 
and X = 0.75 (red). Inset shows snapshots for the FM at X = 1.75 (left) 
and X = 0.75 (right). Snapshots are oriented with the eigenvector of the 
highest order parameter defined as 5. 



allowing particles to order in a nematic-like regime 
(preferred by particle-particle attractive and steric 
packing interactions). This change in the critical value 
of X owing to confinement is analogous to the 



Gibbs-Thomson effect that describes the lowering of the 
liquid-solid phase transition temperature owing to surface 
tension (72). 

In conclusion, particle-particle order is lowered because 
of frustration between wall and nucleosomal shape (bulk 
columnar tendencies that are incommensurate with spher- 
ical boundary alignment). This leads on average to more 
isotropic nucleosomal packaging within the spherical 
shell, in contrast to the nematic, columnar or higher- 
ordered arrangements that are expected for the same 
density and interactions in bulk solutions in the absence 
of confining boundaries (38). Interestingly, similar frustra- 
tion between bulk and confining walls was also shown to 
determine the order of other highly confined nematogens, 
including nematic droplettes, or tactoids (7,73,74). 



CONCLUDING REMARKS 

Showing remarkable agreement, both SAXS experiments 
and coarse-grained simulations of SV40 chromatin point 
to high nucleosome density at the center of the capsid. 
Simulations further indicate a weak hypertrophic layer 
close to the capsid internal wall. In contrast, the organiza- 
tion of the nucleosome-free dsDNA in DNA-VLPs shows 
low nucleic acid density at the capsid center, and high 
DNA concentration close to the viral capsid (Figure lc). 
These results indicate that addition of nucleosomes to 
dsDNA leads to chromatin condensation and allows the 
DNA to extend to the center of the capsid, whereas the 
attractive boundary and tethering between nucleosomes 
favor a hypertrophic layer at the capsid wall. Moreover, 
computer modeling suggests that accounting for nucleo- 
somal packaging under confinement is sufficient to explain 
the minichromosome organization seen by SAXS. 

Although experiments of unconfined NCPs (38), as well 
as our simulations, indicate that nucleosomes in the bulk 
at the densities found inside the virus are orientationally 
ordered, the simulations of confined ellipsoids predict lack 
of particle-to-particle alignment (Figure 4 left snapshot). 
Yet, within the wetting layer, ellipsoids align with the 
boundary (Figure 4 right), reflecting their preferential 
interaction with the wall. Overall, the presence of a 
nanometric confining boundary stabilizes the isotropic- 
like state with respect to the nematic-like regime. This 
disordering tendency is caused by the incommensurate 
wall versus bulk symmetry, as well as by the small 
number of particles that restricts the maximum possible 
correlation length (75). As a result, we predict that the 
minichromosome within the SV40 capsid exists above 
the order-disorder transition (X is higher than X c ), which 
we term a 'molten droplet' state. A similar balance 
between opposing forces is found in confined nematogens 
such as in tactoids, where bulk nematic tendencies 
compete with surface alignment 'anchoring' forces 
(7,73,74). 

Beyond understanding the structural characteristics of 
the packaged SV40 minichromosome, this combined study 
underscores the important effect that strong confinement 
may have on the adaptable and malleable character of 
chromatin ordering. Together with numerous important 
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modifiers of chromatin order that have previously been 
studied (30), it remains to be determined whether such 
'crowding' affects chromatin structure in vivo (76). 
Indeed, it will be interesting to find evidence for such 
nematic 'melting' for other viruses, or for polymers and 
liquid crystals that are subject to nanoconfinement. 

SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online: 
Supplementary Figures 1-9, Supplementary Methods, 
Supplementary Results and Supplementary References 
[77-90]. 
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